Operator space entanglement entropy in transverse Ising chain 
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The efficiency of time dependent density matrix renormalization group metfiods is intrinsicaily 
connected witfi the rate of entanglement growth. We introduce a new measure of entanglement in the 
space of operators and show, for transverse Ising spin 1/2 chain, that the simulation of observables, 
contrary to simulation of typical pure quantum states, is efficient for initial local operators. For 
initial operators with a finite index in Majorana representation, the operator space entanglement 
entropy saturates with time to a level which is calculated analytically, while for initial operators 
with infinite index the growth of operator space entanglement entropy is shown to be logarithmic. 
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I. INTRODUCTION 



The entanglement is an intrinsic property of composite 
quantum systems and represents a cornerstone in quan- 
tum information theory []| . It is important to understand 
the role of quantum entanglement in classical manipula- 
tion of quantum objects, and to quantify the degree of 
entanglement. Although the question of quantification 
is not clearly resolved, the quantum information theory 
offers several measures 0, Q of entanglement. Quan- 
tum information theory also gave a new birth or fresh 
interpretation of a class of methods for numerical simu- 
lation of many-body quantum systems which, due to the 
exponential growth of Hilbert space, cannot be manip- 
ulated using exact diagonalization. The methods orig- 
inally known as density matrix renormalization group 
(DMRG) [4] deploy the fact that many degrees of freedom 
are redundant in quantum state description; the system 
is therefore adequately described by taking into account 
maximally entangled components only. Thus, sufhciently 
slow growth of the entanglement is of crucial importance. 
DMRG enjoyed remarkable success in determining the 
ground state properties of large one-dimensional quan- 
tum models, for which the degree of entanglement scales 
at most logarithmically with size @,ll,llli|9, 10], how- 
ever its time- dependent version (tDMRG) is often 
plagued by abundance of entanglement with time evolu- 
tion [l^ . For efficient classical simulation of many body 
quantum dynamics using tDMRG it is required that the 
computational costs grow polynomially in time meaning 
that the degree of entanglement of any quantum object 
which can be represented as an element of a scalable ten- 
sor product Hilbert space (either a pure state, or a mixed 
state/operator, etc) must grow no faster than logarith- 
mically. It was recently shown that this is generically 
not the case for quantum chaotic Ising spin chain in a 
tilted magnetic field where the entanglement entropies 
grow linearly and hence the computation costs increase 
exponentially in time p^ . 

In this paper we shall consider the model of quantum 
Ising chain in transverse magnetic field which is intc- 
grable and an explicit analytical solution exists. Cal- 
abrese and Cardy [l^ have shown numerically that the 



growth of entanglement entropy is linear for evolution of 
pure initial states which are ground states of quenched 
hamiltonians; see also Ref . . However, from efficiency 
of tDMRG for time evolution of local operators [l^ one 
may conclude that entanglement entropy computed in 
the space of operators grows only logarithmically. Here 
we address this problem theoretically using the idea [I^ 
of re-formulating the Heisenberg evolution in an algebra 
of operators in terms of a Schrodinger evolution gener- 
ated by a different - adjoint hamiltonian acting on the 
Hilbert space of operator algebra. We show that oper- 
ator space entanglement entropy saturates in time for 
initial local operators of finite index (precise definitions 
follow) and explicitly compute the saturation values in 
the critical case. Further, for initial local operators of 
infinite index we give accurate numerical evidence and a 
theoretical hint that the growth is logarithmic (in ther- 
modynamic limit) with prefactor 1/6 in critical, or 1/3 
in non-critical case. 

We note that, to best of our knowledge, the entan- 
glement in operator space is a new concept which has 
not yet been considered theoretically, and it is clearly 
not equivalent to the conventional concept of entangle- 
ment of density operators as discussed in sectllVl Yet, 
it is the one which we expect to be more closely related 
to computational complexity of time-evolution in infinite 
interacting quantum systems. 



II. FERMION REPRESENTATION OF 
DYNAMICS IN OPERATOR SPACE 



The dynamics of a transverse Ising chain of length 2L is 
described in terms of canonical Pauli matrices cr" for sites 
j G '^2L = {^L + 1, . . . , 0, 1, . . . , i} and the hamiltonian 
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with open boundary conditions, which can be diago- 
nalized by means of Jordan- Wigner transformation and 
introduction of Majorana fermion operators [l^, [l7| . 
- ( n,<„ f^lX and r„ = ( n,<„ 'J^H ™iiig the 
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anti-commutation relations {Xi,Xj} = {Yi,Yj} = 2Sij the adjoint reduced density matrix 
and {Xi,Yj} — 0. Heisenberg equations of motion 
dA/dt — i[H, A] for Majorana operators can be written: 



Xn — '2{Yn^i — hYn 



y„ = -2(x„+i-/iX„). (2) 



An operator corresponding to an arbitrary physical ob- 
servable can be written as a superposition of prod- 
ucts of Majorana operators Xj, Yj, namely Pn.n' = 

XlL+iy'-L+i ■ ■■Xl'Y^ with powers n,,n'^ G {0, 1}. 
A set of 4^^ operators {|Pn,n')u} spans an orthonormal 
basis of a Hilbert space, namely the matrix alebra 21 = 
(^2^^x22^^ with an inner product j(A|B)tt = 2-'^^trA^B. 
21 can be conveniently interpreted as a Fock space of 2L 
adjoint fermions (we shall call them a-fermions) with 
pseudo-spin (distinguishing between Majorana Xj and 
Yj operator). An arbitrary operator A is then a-fermion 
state = J2n n' '^n,n'|^n.n')j- A-fermi operators over 
21, Cj|,Cj|, can be introduced by 



Cj| |-Pn,n')jj — nj\XjPn.n')^ 



n',\Y,Pn^n, 



(3) 
(4) 



satisfying canonical anti-commutation relations. 

The index of an operator Pn,n' is defined as /n,n' — 
'^j{nj +nj) and for index-one operators eq. ([2|) rewrites: 

^|^n)„ = 2{cl_,^-hcl)c„^\X„)^ (5) 



-|y„)„ = -2{cl^,. - hclJc^^\Y„ 



dt 



(6) 



which can be interpreted as a Schrodinger equation 
(d/di)|A)u — —iH\A)a for the adjoint hamiltonian 



n = 2i 



(7) 

Since the adjoint time-evolution is a homomorphism the 
Schrodinger equation extends to arbitrary element of the 
operator algebra |A)u G 21. Note that the number of 



-fermions, — ^ 



St 



n,s n,s 



Cn,s is conserved, unlike the 



number of ordinary Majorana fermions. 



III. OPERATOR SPACE ENTANGLEMENT 
ENTROPY 

It is clear that classical simulability of quantum states 
is quantified by the entanglement entropy of half-half (or 
worst case) bipartition of the lattice. However, for simu- 
lability of quantum observables (or density operators of 
mixed states), the decisive quantity is an analog of en- 
tanglement entropy defined for an arbitrary element of 
operator algebra 21 9 \A)a a„,„' |Pn,n')j4, with 



-R(n_i^i,n'_^j,...no,np),(m_L+i,m'_j^j^,...mo,mo) — 

^ ^ 1(n_i^i,...no,ni,...nL),(n'„j^i,...nf|,n'j^,...n^) ^ 

(8) 



namely 



S = -trRlnR. 



(9) 



For a spin 1/2 chain it is perhaps more natural to use 
a set of 4^^ Pauli operators Qs_^ 



-lA-l 



'L ' 



where Sj G {0, x, y, z}, cr" = 1, as a physical basis of oper- 
ator algebra 2t, and define bi-partition and entanglement 
entropy with respect to Qs- However, it is easy to show 
that the result is identical to Q since the transformation 
between the bases {Pn.n'} and {Qs} is a simple permu- 
tation of multiindices (n, n') ^ s (with multiplications 
by ±1), and even though it is non-local it maps first L 
a-fermions to only first L spins and vice versa. 

Let us now try to compute time- dependence of op- 
erator space entanglement entropy S{t) for some sim- 
ple initial operators A. For convenience, we introduce 
staggered a-fermi operators Wj, j G Z4L = {—2L + 
1, . . . , 0, 1, . . . , 2L}, such that 'W2n~i = Cni and W2n = 
Cni- Any operator acting solely in a space of first L 
a-fermions (or first L spins) can be expressed in terms 
of 2L anti-comuting operators wj with j G Z^j^ = 
{-2L + 1, . . . , -1, 0}. We follow Ref.[l3] and express 
2^^ eigenvalues of adjoint reduced density matrix R, as 
Pn = Uj {^Hi + (l-"i)(l-7i)), G {0, 1} where 7^ 
are eigenvalues of time-dependent 2L x 2L correlation 
matrix 

rm„(t) = s(A|w1'„(i)u;„(t)|A)j, m,nGZ-^ (10) 
Then, the entanglement entropy ([5]) simply reads 
S'(i) = ^e(7j), e(a;) = -xlnx- (l-x)ln(l-a;). (11) 

This procedure (see [3| for details) results in an efficient 
numerical method which essentially only requires diago- 
nalization of 2L dimensional matrix F for the solution of 
a quantum problem on 2^^ dimensional Hilbert space. 

The time-dependent a-fermi operators Wm{t) are ob- 
tained from linear Heisenberg type equations Wm = 
—i[wm,'H], namely, W2j = 2{w2j+i — hw2j-i) and 
W2j-i = 2[—'W2j-2 + hw2j)- The solution of such Heisen- 
berg equations, written as Wm = —iJ^n^-mnWn, is ob- 
tained by diagonalizing a.2L x 2L matrix G = V ■ A ■ V"^ 
which yields 
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V, 



71 k 



(12) 



However, in the 'critical case' h — I, the time-evolution 
of Wm{t) can be solved exactly and some analytical so- 
lutions can be given. Namely in such a case the sets of 



Heisenberg eqs. are identical and are solved via discrete 

1^ (m+2L)fc7r] 



sine-transform with V, 



7nk 



4L + 1 



sm 



4L + 1 



4L 

E 

fc=0 



71 k 



(13) 



In the following we shall be interested in the results in 
the thermodynamic limit (TL), L — > oo. The infinite sum 
over k in l|13p is transformed onto an integral which yields 
Wm{t) — X^nez ™ tcmis of Bcsscl functions 

<^ab{x) = Ja-b{x). The Correlation matrix elements are 
therefore (using hb 



wlm) 



rm„(t) =J2^bm{4.t)^bn{it)i{A\nb\A)^. 



(14) 



We also assume that the initial operator A is local, i.e. 
a product of finite number of Pauli matrices cr". This 
implies that: (i) either A has a finite index, i.e. {nb)^ = 
^{A\nb\A)^ = 0, for |5| > bo for some bo £ Z+, or (ii) 
A has an infinite index and {hb)f, = 1 for b < —bo and 
(nb)(t = for & > 6o (such as e.g. A = cr^). Then as 
L — > oo, an arbitrary large fixed finite piece of correlation 
matrix can be asymptotically written as 



Jb^m{^t)Jb-n{^t){nb)i. (15) 



Note that Vmn has effectively finite rank 



4i, 



namely F,, 



for —m,,—n > x. For brevity we 



shall be omitting the argument of Besscl functions al- 
ways equal to x = 4t. 



A. Initial operators of finite index 

First, we focus on the case (i) of finite index initial 
operators A. It was conjectured in [13] that in such 
cases the entanglement entropy in thermodynamic limit 
saturates in time. Using the a-fermion algebra we are 
now able to calculate the exact saturation value since 
RHS in (|15p is a finite sum. Consider F,„„ as a real 
matrix over M°° with canonical basis {\m,),m £ Z~} 
and write a set of non- orthogonal vectors {jV^Q,)}, namely 
{m\ipa) = (— l)"</Q-m(4i) = {ipa\m). Let us write ini- 
tial operator of finite index (K) as A — Oj^ ■ ■ ■ Oj,^ 
where 02j-i = Xj and 02j = Yj. Then we have 



J. 



■JK 



T„rn{t) = (m|V'il)(V'ill?^) + • • • + (m|V'j^)(V'j^|n). (16) 

This means that the rank of Tmn is bounded by K, in 
fact it is K, and its non-trivial eigenspaces are spanned 
by {IV'jJ,! < k < K}. Let {|(/)fc),l < fc < i^} be an 
orthonormalized set obtained from {\ipj^),l < fc < K} 
by a standard Gramm-Schmidt procedure, for which 
the only input is the set of scalar products (V'alV'/s) = 
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FIG. 1: Entanglement entropy for finite-index operators \Xi)^ 
(black) and IXiYi)^ — i|crf)j| (gray) compared to theoretic 
saturation value for t — > oo and h — 1 (thick lines). Three 
different values of magnetic field are considered: h — 1 (solid 
curve), h — 0.5 (dotted), h = 2 (dashed), h — 3 (dash-dotted). 
We set 2L = 200, such that no finite size effect were noticable. 



J2k&z- Jk-aJk-p which Can be in TL evaluated analyt- 
ically for any t in terms of finite sums and approach 
the following long time asymptotics {4'a\'4'p)\t=oo = 
(l/2)(5a/3 — sin[7r(a — /9)/2]/[7r(Q; — /?)]. The non-vanishing 
part of the spectrum {7^} of the correlation matrix (jl6p 
entering eq. pT|) is thus given by the eigenvalues of the 
following K X K matrix 



ffei = {<l>k\ip3i){^n\(l^i) 



{<l>kHj.){i^jM- (17) 



Thus, eq. (|17p is our main result for the case of fi- 
nite index initial operators. For illustration, let us cal- 
culate the asymptotic value S{t — > 00) for the sim- 
plest two cases: (a) A — Oj, e.g. A = ■ ■ ■ a't2<^tiO'o ^ 
and (b) A = OjOj+i, e.g. A = af. In case (a), 
K — 1, the result is 71 — with 7i|t=oo — \ 

which gives the entanglement entropy of S{oo) — In 2. 
In case (b), X = 2, we have 71^2 = ^[(^il^i) + 



+ 4(V',#,+i)^] 
21n(7r'^^72"^=)- 



{ipj+i\i!j+i)± ^J{{ipj\i}j) + (V'j+ilV'i+i))' 

with 7i,2|t=oo = \ ± \ and S'(oo) = 
Both results agree excellently with numerical solutions 
of Heisenberg eqs. for 'Wj{t) shown in Fig.[l] for the case 
h = 1, whereas saturation is observed for any h. 



B. Initial operators of infinite index 

Second, we consider the case (ii) of infinite index ini- 
tial operator A. In thermodynamic limit, local spin 
operators such as af are products of infinite number 
of Majorana operators X„, Yn, in particular |crf)jj = 
I • • • X-iY-iXqYoXi)^ and previous disciussion does not 
apply. As conjectured in time complexity for such 
initial operators only grows polynomially in time which 
corresponds to logarithmic growth of the entanglement 
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FIG. 2: Entanglement entropy for infinite-index operators 
|-F)j| (black) and l^^yL)!] = I'^'i))! (gray) and different mag- 
netic fields (same styling as in Fig. [ij as they appear in the 
legend from top to bottom. Thick dashed lines correspond to 
(1/3) In t and (1/6) Int. 



entropy. Let us define an infinite index operator F = 
• • • X^iY^iXi^Yi^ which corresponds to a half-filled Fermi 
sea \F)^ of a-fermions. Any operator of interest here can 
be written as A = FB where _B is a finite-index opera- 
tor, again can be intepreted as a Fermi sea with a 
finite number of particle and hole excitations. Figure [2] 
shows results for S{t) (fTT|) based on numerical solution 
of Heisenberg equations for Wj (t) for 2L up to 800 such 
that no finite size effects are noticable in the figure. For 
any initial operator of the form A = FB, we observe a 
clean asymptotic logarithmic growth 



S{t) X c\nt + c' , where c : 



1/6 
1/3 



if |/i| = l; 
if\h\^l, 



(18) 



and c' is a constant which, for given h, only depends on 
the choice of finite index operator B. Note an intriguing 
similarity with the size scaling of entanglement entropy 
of the ground state of ([J) [T^ ■ Analytical explanation 
of this interesting phenomenon may be non-trivial, so in 
the following we limit ourselves to the case of critical field 
h = 1 and consider only the simplest initial operator of 
infinite index, namely A = F where the problem can be 
connected to the theory of block Toeplitz determinants. 
In order to avoid negative indices, we re-define the corre- 
lation matrix as rj„„ = r_m_„„, so from eq. follows 



6=0 



m, n = 0, 1, . 



(19) 



It should be noted that correlation matrix can be fac- 
torized T'^„ — J^'iZo^mi'^in as a square of a matrix 



Ho Hi 



V 
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n. 



J21 



■hi+i 
-J21 



(20) 



Following Ref. ^IQ] we express the time-dependent en- 
tanglement entropy (jlip in terms of a formula involving 
Block Toeplitz determinant 



S 



2^ X'(^')[^^''^^*(^^~*)]^^ 



(21) 



where S is a closed curve in complex plane enclosing unit 
disk, avoiding point 1 and interval [—1,0]. Note that 
eigenvalues of infinite dimensional matrix \E' come in pairs 
A, —A with accumulation points ±1. For any e > there 
is only a finite number, Ne{t) ^ t, of eigenvalues of ^ 
which are not in e vicinity of ±1. However, we find nu- 
merically that most of these eigenvalues cluster around 
0, and only ^ In t of them lie outside e vicinity of which 
are the only eigenvalues contributing to entanglement en- 
tropy result 

At present state of the theory of block Toeplitz 
determinants - in connection to the theory of inte- 
grable Fredholm operators and the Riemann-Hilbert 
problem [l9| - the formula ([2T|) can be explicitly eval- 
uated (see e.g. Ref. [20j ) provided the matrix sym- 
bol ^{z) = Al — X^fcez^fc^'^j admits explicit Wiener- 
Hopf factorizatons $(z) = C/+(z)C/-(z) = V~{z)V+{z) 
where the matrix functions U^{z),V^{z) are analytic 
inside(-f )/outsize(— ) the unit circle. Even though the 
matrix symbol has an appealing explicit form 



$(z) = 



A - // + gg fg/ z- gf 
zgf - fg + ff - gg 



(22) 



(— l)"J„i-„. Note that is a real symmetric 



where / = f{z)J = f{z ^),g = g{z),g = g{z ^) and 
f{z) = cosh{2t^/z) and g{z) = smh{2ty/z) / ^/z are en- 
tire analytic functions, its Wiener- Hopf factorizaton is at 
present unknown and poses a challenging problem. 



IV. CONCLUSIONS 

We have studied complexity of time evolution of ini- 
tial local operators under dynamics given by the trans- 
verse Ising chain. Such complexity can be characterized 
in terms of entanglement entropy of operators treated 
as elements of a product Hilbert space corresponding to 
a bi-partition of a chain and is directly related to time 
efficiency of simulation methods such as tDMRG. Note 
that operator space entanglement entropy, of say a den- 
sity operator, is not simply related to a traditional notion 
of entanglement of the corresponding mixed state. For 
example, consider a macroscopic convex combination of 
2^ product states. This corresponds to a non-entangled 



5 



mixed state but has a macroscopic operator space entan- 
glement entropy ^ L and hence it is difScult to simulate 
classically. Thus it seems that the traditional concept of 
state entanglement is not sufficient to characterize clas- 
sical complexity of quantum operators. In this paper we 
have shown, in parts analytically and numerically, that 
operator space entanglement entropy of transverse Ising 
model grows at most logarithmically for initial operators 



which are local products of Pauli matrices. This result 
has to be contrasted with a a linear growth of entan- 
glement entropy for time evolution of pure states [3]. 
Explanation of deeper physical reasons for this dramatic 
effect is needed. 
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